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Abstract. In this paper we model and analyse the B — V/V ~ K radial color gradients observed in the Cartwheel 
ring galaxy. Along with the color-color diagrams, we use the Qbvk combined color indices, which minimise the 
uncertainties in the observed B — V and V — K colors introduced by dust extinction. To model the optical and 
near-infrared color properties of the Cartwheel galaxy, we assume that an intruder galaxy generates an expanding 
ring density wave in the Cartwheel's disk, which in its turn triggers massive star formation along the wave's 
perimeter according to the Schmidt law. We use the population synthesis to calculate the color properties of 
stellar populations formed in the expanding density wave. The results of color modelling suggest that the pre- 
collision Cartwheel was a late-type spiral, embedded in an extensive gaseous disk of sub-critical surface density. 
The properties of the old stellar disk are typical for the late-type Freeman disks, with the central surface brightness 
in V-band and the scale length being fi v — 21.0 mag arcsec -2 and Ro = 2 kpc respectively. The pre-collision 
gaseous disk has a metallicity gradient ranging from z — zq/5 at the outer regions to z — zq in the central regions. 
At present, the wave of star formation has passed the initial extent of the pre-collision, old stellar disk and is 
currently moving in the predominantly gaseous, low-metallicity disk at the radius of 16 kpc. Neither young stellar 
populations formed in an expanding density wave, nor their mixture with the old, pre-collision stellar populations 
can reproduce the B — V and V — K colors of the Cartwheel's nucleus+inner ring. We find that an additional 
10-Myr-old burst of star formation in the nuclear regions, along with the visual extinction of Ay = l m .3, might 
be responsible for the peculiar colors of the Cartwheel's nucleus. 
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1. Introduction 

Ring galaxies are believed to be the result of a head-on 
galaxy-galaxy collision. Such a near-central collision gen- 
erates an expanding ring density wave in the disk of a 
larger galaxy. It is likely that star formation will be in- 
duced along the wave's perimeter when the gas density 
exceeds a threshold (Appleton & Struck-Marcell 1987). 
As such a wave of star formation advances radially from 
the nucleus, it leaves behind evolved stellar populations, 
with the youngest stars located at the current position of 
the wave. It is expected that this might result in the ra- 
dial color distribution of stars in the galactic disk, with 
the inner regions being redder than the outer parts of the 
disk. Indeed, o bserva tions of the Cartwheel ring galaxy 
(Marcum et al. 1992j ) reveal the presence of the optical 
and near-infrared B — V/V — K radial color gradients in 



the galactic disk. Recent observations of a samp le of north- 
ern ring galaxies by Appleton & Marston ( |l997|) show that 
most of them exhibit radial color gradients similar to the 
pattern observed in the Cartwheel. 

Numerical modelling of the Cartwheel's radial color 
gradients by Korchagin et al. ( |200l| ) shows that the young 
stellar populations formed in the expanding density wave 
exhibit regular reddening of colors towards the nucleus 
only for the sub-solar metallicities of the star-forming gas 
(z < zq/2.5). Even in the most favorable case of z = Zq/5, 
the model colors are much bluer as compared to those ob- 
served in the Cartwheel. Korchagin et al. (2001) argue 
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that the pre-collision disk of old stars is needed to recon- 
cile the model and the Cartwheel's colors. However, their 
conclusions are based on the direct quantitative compari- 
son of the model and the Cartwheel's B — V and V — K 
colors, which is complicated due to uncertainties in the 
amount of internal extinction in the Cartwheel's disk. 
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In this paper we re-address the question of theoreti- 
cal modelling of the B — V/V — K radial color gradients 
observed in the Cartwheel galaxy. We attempt to exclude 
the complicated influence of dust extinction on the the- 
oretical interpretation of the Cartwheel's B — V/V — K 
radial color gradients. To do this, we use insensitive to 
dust extinction Qbvk combined color indices. Such color 
indices were successfully applied for structure modelling 
of dusty galaxies by Zasov & Moiseev (1998) and Bizyaev 
et al. ( pOOOt pOOll ). 

In Sect. 2 we analyse the Qbvk radial distribution ob- 
served in the Cartwheel's disk and consider the factors 
defining the value of Qbvk- In Sect. 3 we formulate the 
adopted model of star formation in the Cartwheel galaxy. 
In Sect. 4 we compare the model B-V/V-K radial color 
gradients and Qbvk index profiles, obtained in the frame- 
work of a first major episode of star formation in the 
Cartwheel's disk, with the observed radial profiles. We 
consider a possibility that the Cartwheel has an old stellar 
disk typical for the late- type spirals. In Sect. 5 we attempt 
to model the peculiar colors of the Cartwheel's nucleus 
and the inner ring. The main results are summarized in 
Sect. 6. 
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2. Infrared and optical B — V/V — K radial color 
gradients and Qbvk combined color indices of 
the Cartwheel galaxy 

The B — V/V — K radial color gradients of the Cartwheel 
galaxy are shown in Fig. 1, which is a reproduction 
of Fig. 3 of Marcum et al. (1992). The annular points 



are numbered according to radius, beginning with the 
inner ring (point I) and ending with the outer ring 
(points VIII and IX). These colors have only been cor- 
rected for Galactic reddening using A-q = m .22 and the 



standard ISM extinction curve (Marcum et al. 1992). 

Previous modelling of the Cartwheel's radial color gra- 
dient was seriously complicated by the uncertainties in the 
amount of internal extinction in the galactic disk. The in- 
ternal extinction is only measured towards two bright HII 
regions in the Cartwheel's outer ring (Ay — l m .94) and 
is not applicable to the whole galactic disk. The value 
of internal extinction interior to the outer ring of Ay = 
m .3 has been estimated from the available HI + H2 data 
(Korchagin et al. |2001 ). Clearly, correction for the internal 
extinction would greatly distort the observed Cartwheel's 
B-V/V-K radial color pattern, as shown by correction vec- 
tors in Fig. 1, and complicate the task of color interpreta- 
tion. 

To avoid the complicated influence of dust on the 
Cartwheel's radial color gradient, we use insensitive to 
dust extinction Qbvk combined color index defined as: 



hvK = (B-V)-^^(V-K), 



(1) 



where (B — V) and (V — K) are colors, and E(B — V) and 
E(V — K) are color excesses. 



Fig. 1. Radial B — V/V — K color gradients of the 
Cartwheel gala xy. T his is a reproduction of Fig. 3 of 
Marcum et al. ( |l992|) . 



For dust properties typical for the Galaxy (Ry = 3.1) 
the ratio of color excesses E(B — V)/E(V — K) is equal 
to 0.365 (Cardelli et al. |1989| ). This value is valid for 
a "dust shield" approximation, i.e. when all dust is lo- 
cated between an observer and a galaxy. For a homo- 
geneous mixture of dust and stars the ratio of color ex- 
cesses E(B - V)/E(V - K) is 10% less than for a "dust 
shield" approximation (Zasov & Moiseev 1998| ) and equal 
to 0.329. As to the Cartwheel galaxy, the internal extinc- 
tion dominates at the outer ring and is comparable to 
the Galactic reddening in the inner part (Korchagin et al. 
pOOlt ). Hence, we choose a "homogeneous mixture " ap- 
proximation rather than a " dust shield" for modelling the 
Cartwheel's color properties. 

The following factors define the value of Qbvk : 
a) contribution of young, blue stars to the integrated 
flux, i.e. the history of star formation in a certain region 
of a galaxy. Figures 2a and 2b show the Qbvk combined 
color indices of a single starburst superimposed on the old 
stellar population of 10 Gyr. Calculations are performed 



using the program of Worthey (1994). Metallicities of both 
the young and old stellar populations are equal to Zq/5 
([Fe/H]=-0.7). Figure 2a illustrates the influence of dif- 
ferent relative contributions of young stars on the Qbvk 
indices for a single 50-Myr-old starburst. Figure 2b shows 
the dependence of the Qbvk indices on the age of a star- 
burst. Here, both the young and old stellar populations 
contribute equally to the total stellar mass. As it can be 
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Fig. 2. The Qbvk indices of a single starburst superim- 
posed on the old stellar population of 10 Gyr. a) The 
influence of different relative contributions of young stars 
on the Qbvk indices for a single starburst of 50 Myr old. 
b) Dependence of the Qbvk indices on the age of a star- 
burst, c) The Qbvk indices as a function of mctallicities 
of old and young stellar populations (equal for both pop- 
ulations) . 



seen, Qbvk indices are sensitive to the age of a starburst 
and to the relative contribution of young stars. 

Thus, a young burst of star formation should essen- 
tially influence the value of Qbvk- Relative increase of 
young stars in a certain region of a galaxy makes its colors 
bluer, thus, lowering the value of Qbvk (Fig- 2a). On the 
other hand, the red supergiant flash (Chariot & Bruzual 



1991) makes the luminosity of a young stellar cluster in 



the if-band to peak at 10 Myr after the burst, which no- 
ticeably lowers the value of Qbvk (Fig. 2b). 

b) heavy element abundances. Figure 2c gives Qbvk 
indices as a function of metallicities of the old and young 
stellar populations (equal for both populations). The age 
of a single starburst is equal to 50 Myr. Metallicity gra- 
dient alone can essentially influence the value of Qbvk, 
while other starburst parameters remain unchanged. 

c) dust properties. Cardelli et al. ( |l989|) found that for 
the dust properties different from those of the Galaxy, i.e. 
for i?y other than 3.1, the ratio of color excesses E(B — 
V)/E(V — K) changes. However, the Cartwheel's Qbvk 



distribution remains qualitatively unchanged for a fairly 
wide range of Ry. Thus, the value of Ry — 3.1 is adopted 
for the rest of the paper. 

Figure 3 shows the Cartwheel's Qbvk distribution as 
a function of photometric annulus. Combined color in- 
dices are calculated using observational data of Marcum 
et al. (1992). The numbering of photometric annuli is 



in accordance with Fig. 3 of Marcum et al. (1992). The 



Cartwheel's Qbvk radial distribution remarkably sepa- 
rates into three sub-profiles, which correspond to three 
structurally different parts of the Cartwheel: annuli and 
1 correspond to the nucleus and the inner ring, annuli 2-7 
are the inter-ring region, and annuli 8 and 9 correspond to 
the outer ring. Combined color indices are insensitive to 
dust extinction and can be directly understood in terms of 
stellar populations. This natural separation of Qbvk in- 
dices implies that the nuclcus+inncr ring, inter-ring region 
and the outer ring of the Cartwheel consist of different 
stellar populations. 

Nucleus and the inner ring. The old stellar content 
of the nucleus is confirmed by the detection of lat e-typ e 
stellar absorption spectrum (Fosbury & Hawarden 1977). 
However, recent mid-infrared and Ha observations indi- 
cate that the old stellar populations may not totally dom- 
inate the spectrum. ISOCAM images (Charmandaris et 
al. 1999] ) show significant 7 \im and 15 fj,m emission in the 
nucleus and the inner ring, which implies a low-level star 
formation activity. Detection of weak Ha emission from 
the nucleus and inner ring also indicates that there might 
be some low-level star formation in the Cartwheel's central 
regions (Amram et al. 1998). 

The inter-ring region. As Ha observations indicate, 
it is almost totally devoid of any star formation activ- 
ity (Amram et al. 1998) and is probably populated by 



density wave born, post-collision stars with an admixture 
of old, pre-collision stars as implied by progressively red- 
dened B—V/V—K colors towards the Cartwheel's nucleus. 
The possibility that dust extinction is responsible for this 
reddening of colors will be tested in numerical simulations 
in Sect. |. 

The outer ring. Models of ring galaxy formation 
predict that the Cartwheel's outer ring is the current 
location of an outwardly propagating star formation 
wave (Ap pleton & Struck-Marcell |l987j , Struck-Marcell & 
Appleton 1987 ). Hence, a large gradient in the Qbvk dis- 
tribution seen across the outer ring in Fig. 3 is indicative 
of rapid changes in the photometric properties of young 
stellar populations. Probably, the outer edge of the ring 
is so young that it has not yet produced any supergiants, 
whereas the inner edge is in the red-supergiant-flash stage 
(Marcum et al. 1992). Existence of the red supergiant flash 
is confirmed in observations of the young stellar cluster 



NGC 2004 in the LMC (Bica & Alloin 1987) 
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Fig. 3. The Cartwheel's Qbvk distribution as a function 
of photometric annulus. 



where t,„ is the life-time of a star of mass m, and a is the 
slope of the IMF. 

We consider a circular density enhancement propa- 
gating from the center of a gaseous disk to the present 
location of the Cartwheel's outer ring at 16 kpc (Ha = 
100 km s _1 Mpc -1 ). The computational area is divided 



into 10 annuli as in Plate 1 of Marcum et al. ( 1992 ). 
Annulus 1 is two kiloparsec wide and centered at 3.1 kpc, 
thus covering the current location of the inner ring. The 
centroids of annuli 2-9 are separated by 1.25 kpc, with 
the 2-nd annulus centered at 6.8 kpc. The 0-th annulus 
represents the nucleus. The mass of stars formed in each 
annulus is distributed into individual stellar masses using 
the Salpeter IMF with a = 1.35 and the stellar mass in- 
terval of mi=0.1 Mq and m2 = 100 M Q . The luminosity 
of each annulus at a given band A after time t can be 
computed as: 



L A(t) = X! X! l A{m, T m )N(m, T m , t), 



(5) 



3. Theoretical modelling of the Cartwheel's color 
properties: orbit crowded density wave plus a 
Schmidt law of star formation 

To model the optical and near-infrared color properties of 
the Cartwheel galaxy, we assume that the head-on galaxy- 
galaxy collision generates an outwardly propagating ring 
density wave in the disk of the Cartwheel. Additionally, we 
assume that the collision is central or near-central so that 
we can neglect the azimuthal variations in the strength of 
the density wave. Massive star formation in the Cartwheel 
galaxy is restricted to the outer ring, which is the current 
location of the density wave. Almost no sigh of star forma- 
tion activity is found just behind and ahead of the outer 
ring (Higdon 1995 ). 

To imitate the propagating star formation in the 
Cartwheel's disk, we assume that the star formation rate 
is proportional to a ring density enhancement expanding 
radially in a gaseous disk with velocity v. The density en- 
hancement is given by a Gaussian function: 



C(r — vt) = A(r) exp 



(r - vt) 2 



L 2 



(2) 



where A(r) is the amplitude of the density enhancement 
and L is a half- width of the density enhancement. The 
balance of surface density of stars (S s ) at a distance r 
from the center and time t can be then written as follows 



(Korchagin et al. 2001): 



^ = -D + aC n {r-vt), 
at 



(3) 



where the first term in the right-hand side of Eq. is the 
rate of death of stars: 



D=- J2 ml ^ aCn ( r - v ^- T ^^ ( 4 ) 

mi 



,1-a 



where N(m,T m ,t) is the number of stars in a given an- 
nulus at time t with mass m and age r m , and lA(m,T m ) 
is the luminosity of a star in the band A. Stellar lumi- 
nosities lA(m,T m ) are obtained using stellar evolutionary 
(Schallcr et al. 1992| ) and atmospheric (Kurucz 1992) mod- 
els. The population synthesis code used in this work is 
explained in detail in Mayya (1995, |1997 ). This code syn- 
thesizes a number of observable quantities in the optical 
and near-infrared parts of the spectrum, which are suit- 
able for comparison with the observed properties of giant 
star-forming complexes. Since the Cartwheel is a recent 
phenomenon with high rates of star formation, this code 
is especially suitable for our purposes. The results of this 
code are compared with those of other existing codes by 
Chariot ( JT99^ ). 

Classic Schmidt law of star formatio n (n = 1.5) is as- 



sumed throughout the paper (Kennicutt 198S ). The half- 
width of a density enhancement L is chosen so as to re- 
produce the observed Ha surface brightness profile in the 
Cartwheel's outer ring (Higdon |i"995l ). This gives the value 
of L = 1 kpc. Velocity of a propagating density enhance- 
ment v = 55 km s _1 is assumed to coincide with the ex- 
pansion velocity of the outer ring v = 53 ± 9 km s , the 
value derived by Higdon ( 1996 ) from HI kinematics of the 
outer ring. 

With the values of L and v being fixed, the coefficient a 
and the amplitude of density enhancement A(r) constitute 
the set of parameters, which determine the number of stars 
N(m, T m , t) formed by a density wave in each photometric 
annulus (see eq. (||)). The birthrate of stars a relates the 
star formation rate observed in the Cartwheel galaxy with 
the assumed Schmidt law of star formation (see eq. (|J)). 
We estimate its value for the Cartwheel's outer ring using 
the measured rate of star formation 67 Mq yr _1 (Higdon 



|l995 ), t he me an surface density of gas S gas w 10 M Q pc~ 2 
(Higdon 1996| ) and the outer ring's area of 250 kpc 2 . This 
gives the dimensionless value of a = 0.0025 when masses, 
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time scales, and distance scales are expressed in units 
10 7 Mq, 10 6 yrs, and 1 kpc respectively. 

We assume that the amplitude of density enhancement 
A(r) is twice as big as the unperturbed, pre-collision gas 
surface density, which is typical for the spirals and is con- 
firmed in numerical simulations by Appleton & Struck- 
Marcell ( 1987 ) . To constrain the pre-collision radial profile 
of gas, we fit the model i?-band and Ha surface bright- 
ness profiles of stellar populations formed in the expand- 
ing density wave to the profiles observed in the Cartwheel. 
We do this by varying the radial profile of gas in the pre- 
collision galactic disk. 

4. Model B-V/V-K radial color gradients and 
Qbvk indices of the Cartwheel galaxy. 

4.1. Density wave in a purely gaseous disk. 

As mentioned in the introduction, Korchagin et al. ( [2001 ) 



found that the density wave propagating in a purely 
gaseous disk of uniform metallicity cannot reproduce the 
Cartwheel's B — V/V — K radial color gradient. The model 
colors are much bluer as compare d to those observed in 
the Cartwheel. Korchagin et al. ( 2001 ) argued that the 
pre-collision disk of old stars was needed to reconcile the 
model and the Cartwheel's colors. We mainly confirm their 
conclusion. However, use of the extinction-free Qbvk in- 
dices allows us to notice some intere sting features that 
have been missed by Korc hagin et al. ( 2001 ). 

As in Korchagin et al. ( 2001 ), we consider the expand- 
ing density wave in a purely gaseous disk with five differ- 
ent metallicities: Z0/2O, z©/5, Zq/2.5, zq, and 2z©. Our 
computations show that the model B — V and V — K col- 
ors are bluer as compared to the Cartwheel's colors for 
each metallicity of the star-forming gas. If dust extinc- 
tion were solely responsible for this difference between the 
model and observed colors, then the model Qbvk indices 
would lie within the observed limits. Indeed, use of the 
Qbvk indices minimises the uncertainties in observed col- 
ors introduced by dust extinction. However, we find that 
the model Qbvk indices lie beyond the observed limits 
for most of the photometric annuli, thus indicating that 
the reddening of the observed colors with respect to the 
model colors of those annuli cannot be attributed to the 
extinction gradient in the Cartwheel's disk. 

In case of Zq/2.5 and z©/5, however, the Qbvk in- 
dices of the Cartwheel's outer region (annuli VI-IX) lie 
within the observed limits, thus indicating that dust ex- 
tinction could be responsible for the reddening of the ob- 
served colors in this region. Here, we note that Korchagin 
et al. (2001) used a uniform metallicity of Zq/20 to model 
the Cartwheel's radial color gradient. While this choice is 
certainly just ified f or the Cartwheel's outer ring (Fosbury 
& Hawarden 1977 ), it is hard to believe that such a low 
metallicity is present in the Cartwheel's central regions. 
Due to very weak line emission, there have been no mea- 
surements of heavy element abundances interior to the 
outer ring. However, it is known that most galactic disks 



have a metallicity gradient, with metallicities in the nu- 
cleus up to ~ 10 times higher than in the outer parts of the 
disk (Smartt & Rollcston 1997). In this paper we assume 
that metallicity gradient is present in the pre-collision 
gaseous disk of the Cartwheel galaxy and choose its value 
so as to provide a better fit of the model Qbvk indices 
with those observed in the Cartwheel's disk. The best fit 
is found for the metallicity gradient spanning the range 
from z = Z0/3.75 at annuli IX, VIII, and VII, z = z Q /2.5 
at annuli VI and V, to z = zq at annuli IV-nucleus. 

The open circles in Fig. 4a show the model B—V/V—K 
colors of stellar populations formed in the density wave 
propagating in a purely gaseous disk with the metallicity 
gradient given above. The open circles in Fig. 4b illustrate 
the corresponding model Qbvk indices. The outer three 
photometric annuli (9,8, and 7) and the nucleus are iden- 
tified both for the model colors and the colors observed in 
the Cartwheel galaxy, the latter being shown by the filled 
triangles with error bars. The model B — V/V — K colors 
do not exhibit a regular reddening towards the nucleus 
and are bluer as compared to the Cartwheel's colors. On 
the same time, the model and the Cartwheel's Qbvk in- 
dices show a good correspondence for the outer ring and 
two inner annuli adjacent to the outer ring, implying that 
the difference of the corresponding B — V/V — K colors in 
Fig. 4a might be related to the internal extinction gradient 
existing in the Cartwheel's disk. 

Now we consider an extreme case and assume that the 
difference of the model and observed B — V/V—K colors in 
Fig. 4a is solely due to internal extinction gradient in the 
Cartwheel's disk. We warn, however, that this assumption 
is applicable only to those annuli, which Qbvk indices lie 
within the observed limits, as in this case their colors are 
reddened along the correction vector. The values of color 
excess E(B — V), implied by this difference, are listed in 
Col. 2 of Table 1. The surface densities of gas S^ c I t +H2 in 
Col. 3 are obtained from Col. 2 using the standard mass- 
to-dust ratio and Ry = 3.1 (Bohlin et al. 1978j ). Column 4 
in Table 1 gives the observed surface densities of atomic 
and molecular hydrogen S^^jj^ in each annulus (Higdon 



1996, Horellou et al. 1998) 



Table 1. Extinction-estimated and observed surface den- 
sities of HI and H2 in the Cartwheel 



annulus 


E{B - V) 
(mag) 


Mq PC" 2 


Mq PC" 2 


VI 


ra .13 


5.7 


< 4.5 


VII 


m .07 


3.1 


< 4.5 


VIII 


m .14 


6.1 


9.3-16 


IX 


m .25 


10.9 


9.3-16 



It can be noticed that the extinction-estimated gas 
densities in Col. 3 are below the detection upper limit for 
the outer ring+annulus VII. Hence, the internal extinc- 
tion gradient can indeed be responsible for the difference 
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Fig. 4. (a) The B — V and V — K color-color diagram and 
(b) the Qbvk combined index profiles. The open circles 
show the model profiles for the density wave of amplitude 
A(r) = 70 exp(— r/12) A/ Q pc -2 and velocity v = SSkms" 1 
expanding radially in a purely gaseous disk. The metallic- 
ity gradient of the gaseous disk ranges from z = z Q at 
the nucleus +annuli I-IV, z = z Q /2.5 at annuli V-VI, to 
z = zq/3.75 at annuli VII- IX. The filled triangles with er- 
ror bars show the radial color gradients and index profiles 
observed in the Cartwheel galaxy. 



between the observed and model B — V/V — K colors of 
this region in Fig. 4a. 

Very large values of Ha equivalent widths were found 
by Higdon (1995) in the Cartwheel's outer ring, espe- 
cially in the ring's southern quadrant, which peaks at 
EWn a = 1250 A and exceeds 400 A everywhere. The 
Ha equivalent width (EWn a ) is sensitive to the ratio of 
ionizing, massive (> 10 M Q ) stars to lower mass red giant 
stars and can be used to measure the relative strengths of 
these two populations (Kennicutt 1983). Table 2 presents 
EW}i a of stellar populations born in the expanding ring 
density wave. EWua are computed at the time when the 
wave reaches the present position of the Cartwheel's outer 
ring at 16 kpc. All parameters of Fig. 4 are retained. As 
it is seen in Table 2, the averaged over annuli VIII and IX 
model value of 335 A is close to the azimuthally averaged 
value of 360 A measured in the Cartwheel's outer ring by 
Higdon (1995). The large model and observed EW^ a in 



the outer ring indicate that the Cartwheel is experiencing 
the first major episode of star formation at this radius. 

Table 2. Model Ha equivalent widths 



annulus 


VI 


VII 


VIII 


IX 


EWu a (A) 


0.03 


0.3 


53 


620 



Korchagin et al. (2001) argued that the colors of the 
Cartwheel galaxy cannot be reproduced by a density wave 
propagating in a purely gaseous disk irrespective of its 
metallicity. While this is certainly true for the Cartwheel's 
inner region, we find that the observed B — V/V — K colors 
and Qbvk indices of the outer ring+annulus VII can be 
well reproduced by a density wave propagating in a purely 
gaseous disk of z — z Q /3.75. 

4.2. The pre-collision old stellar disk. 

To reconcile the model and observed B — V/V — K colors 
and Qbvk indices of the Cartwheel's inner region in Fig. 4 
(annuli VI- nucleus), we follow the prescriptions given in 
Korchagin et al. (2001) and consider a possibility that the 
pre-collision Cartwheel had an old stellar disk. We assume 
that the properties of the old stellar disk are typical for 
the late-type Freeman disks (fly — 21.0 mag arcsec -2 and 
Ro = 2 kpc) rather than for low surface brightness galax- 
ies (/Xy = 23.5 mag arcsec -2 and Rq = 5 kpc) as adopted 
by Korchagin et al. (|2001 ). Metallicity of the pre-collision 
gaseous disk is assumed to span the range from z = zq/5 
at six outermost annuli to z = zq at the inner annuli and 
the nucleus. Galactic collision generates an expanding den- 
sity wave, which adds the newly born stellar populations 
to the old stellar populations existed before the collision. 

Figures 5a and 5b show the resulting model Ha and R- 
band surface brightness profiles. Interior to the outer ring, 
the model i?-band profile shows a good correspondence 
to the observed profile, as illustrated in Fig. 5b by the 
open cir cles a nd filled triangles respectively. Measured by 
Higdon ( 1995 ) _R-band surface brightness of the inner ring 
is roughly twice the peak value of the outer ring and this 
feature is well reproduced in Fig. 5b. 

Table 3. Extinction-estimated and observed surface den- 
sities of HI and H 2 in the Cartwheel 



annulus 


E{B - V) 
(mag) 


yiext 
Zj HI + H 2 

M Q pc" 2 


v~»obs 
^HI+Hs 

M pc" 2 


II 


m .07 


3.1 


< 4.5 


III 


m .07 


3.1 


< 4.5 


IV 


m .07 


3.1 


< 4.5 


V 


m .07 


3.1 


< 4.5 


VI 


m .07 


3.1 


< 4.5 


VII 


m .06 


2.6 


< 4.5 


VIII 


m .15 


6.6 


9.3-16 


IX 


m .26 


11.4 


9.3-16 
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Fig. 5. The radial surface brightness profiles in (a) Ha and (b) i?-band, (c) the B — V/V — K color-color diagram, 
and (d) the Qbvk radial profiles. The open circles show the model profiles for the density wave of amplitude A(r) = 
50 exp(— r/15) M pc~ 2 and velocity v — 55 km s" 1 expanding radially in the galactic disk. The central surface 
brightness in the F-band of the old stellar component of the galactic disk is /z v = 21.0 mag arcsec -2 and the scale 
length is i? = 2 kpc. The gaseous component has the metallicity gradient ranging from z — zq at the nucleus +annuli I- 
III to z = Zq/5 at annuli IV-IX. The filled triangles with error bars show (a,b) the measured radial surface brightness 
profiles, (c) the radial color gradients and (d) index profiles of the Cartwheel galaxy. 



The model radial B — V/V — K colors, shown by 
the open circles in Fig. 5c, exhibit regular reddening to- 
wards the nucleus reminiscent of that observed in the 
Cartwheel's disk. The sequence of Qbvk indices observed 
in the Cartwheel's outer ring and the inter-ring region is 



well reproduced, as shown by the open circles in Fig. 5d. 
This implies that the difference of the model and observed 
colors for each photometric annulus can be understood 
in terms of internal extinction. This is certainly not true 
for the Cartwheel's inner ring and the nucleus, as their 
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B — V/V — K colors are not reddened along the correction 
vector with respect to the model colors. This results in a 
substantial disagreement of the corresponding Qbvk in- 
dices. As in Sect. 4.1, we consider an extreme case and 
assume that the difference of the model and observed 
B — V/V — K colors of the outer ring and the inter-ring re- 
gion in Fig. 5c is solely due to internal extinction gradient 
in the Cartwheel's disk. The values of E(B — V) implied 
by this difference are listed in Col. 2 of Table 3. 

It can be noticed that the extinction-estimated gas 
densities in Col. 3 are below the detection upper limit 
in Col. 4 for all annuli. Hence, the internal extinction gra- 
dient can be responsible for the difference between the 
observed and model B — V/V — K colors of the inter- ring 
region and the outer ring in Fig. 5c. 

Comparison of Figs. 4a and 5c shows that the model 
B — V/V — K colors of the outer ring+annulus VII in 
Fig. 5c are not affected noticeably by the old stellar pop- 
ulations. Indeed, Fig. 6 illustrates the ratio of fluxes of 
the old stellar populations to the young density-wave- 
born stellar populations in V-band (Fy[o/dw]) and K- 
band (Fk[o / dw]) as a function of galactic radius. It is seen 
that the old stellar populations dominate at the smaller 
radii (R < 3 kpc) in both bands. At the intermediate 
radii (3 < R < 12 kpc), input from the young stellar 
populations to the total flux in F-band is ten times (on 
average) the input from the old stellar populations, while 
both stellar populations contribute comparable fluxes in 
if-band. The young stellar populations totally dominate 
at the larger radii (R > 12 kpc ) in both bands. This in- 
dicates that the colors of the Cartwheel's outer ring are 
determined by the young density-wave-born stellar popu- 
lations and can be reproduced by a density wave prop- 



Distance (kpc) 

5 10 



agating in a purely gaseous disk, as argued in Sect. 4.1 



We note, however, that the interpretation of colors of the 
Cartwheel's inner region do require the presence of the 
old pre-collision stellar disk, w hich is in agreement with 
conclusions of Korchagin et al. ( [2001 ). 

Thus, theoretical modelling of the color gradients ob- 
served in the Cartwheel galaxy suggests that the pre- 
collision Cartwheel was a late- type spiral (R < 12 kpc) 
embedded into an extensive HI disk . The subcritical gas 
disk of S H i < 11 M Q pc~ 2 (Higdon |1996|) at larger radii 
(i? > 12 kpc) prevented robust massive star formation 
(MSF) and chemical enrichment over the disk's lifetime. 
The passage of a companion galaxy initiated an outwardly 
propagating ring density wave, which triggered high rates 
of MSF along it's perimeter. The density wave has passed 
the original extend of the pre-collision stellar disk and 
is currently moving in the predominantly gaseous, low- 
metallicity disk at the radius of the Cartwheel's outer ring 
(R w 16 kpc). 

Another result of color modelling is a possible presence 
of a large internal extinction gradient across the outer ring 
of the Cartwheel (annuli VIII and IX). Indeed, if the dif- 
ference between the model and observed B — V/V — K 
colors of the Cartwheel's outer ring in Figs. 4a and 5c can 
be attributed to the internal extinction as shown above, 



model 'nucleus' 

model 'inner ring' 




0.001 [ — i — i — i — i — | — i — i — i — h 
model 'nucleus' 

model 'inner ring' 




5 10 
Distance (kpc) 



Fig. 6. The ratio of fluxes in V-band (the upper frame) 
and if-band (the lower frame) of the old, pre-collision stel- 
lar populations to the young stellar populations formed in 
the density wave. All parameters of Fig. 5 are retained. 
The numbering of annuli is in accordance with Fig. 1. 



then E(B — V) in the outer edge of the outer ring (annulus 
IX) is about two times higher as compared to E(B — V) 
in the inner edge of the outer ring (annulus VIII). Models 
of ring galaxies predict that the Cartwheel's outer ring is 
the current location of a star formation wave (Appleton 
& Struck-Marcell |1987| , Struck-Marcell & Appleton |l~987| ) . 
The star formation wave is expected to compress the gas 
ahead by a coherent action of shock waves, expanding su- 
perbubbles, ultraviolet radiation from massive stars, ets., 
thus increasing the value of extinction at the leading edge 
of the wave. Observationally proved presence of such an 
extinction gradient across the outer ring might be a pow- 
erful confirmation of propagating nature of star formation 
in ring galaxies. More sophisticated models of star forma- 
tion in the ring galaxies are needed to investigate these 
aspects. 



5. The inner ring and the nucleus 

In the cloud-flu id m odels of ring galaxies, Appleton & 
Struck-Marcell ( 1987 ) predicted that the infall of gas be- 
hind the ring should lead to a strong compression of gas in 
the central parts. In case of the Cartwheel galaxy, a large- 
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V-K 

2 



Cartwheel's 
nucleus 




model 'nucleus' 




3 4 5 6 7 
Photometric annulus 



Fig. 7. (a) The B — V/V — K color-color diagram, and 
(b) the Qbvk radial profiles. All parameters are the same 
as in Fig. 5, but with an additional 10-Myr-old burst of 
star formation in the Cartwheel's nuclear region. 



scale radial infall of HI was reported by Higdon ( 1996| ). 
Recent observation in the mid-infrared (Charmandaris et 
al. 1999) and in Ha line (Amram et al. 1998) have shown 
that the Cartwheel's inner ring and the nucleus might have 
a low level star formation activity. 

In this section, we perform numerical modelling of the 
B — V/V — K colors and Qbvk indices measured in the 
Cartwheel's inner ring and the nucleus (Marcum et al. 
1992| ). To do this, we assume that the star formation takes 



place not only along the perimeter of the expanding den- 
sity wave, but also in the central region of the Cartwheel. 
We set off an outwardly propagating density enhancement 
of amplitude A(r) = 50 exp(— r/15) Mq pc~ 2 and veloc- 
ity v — 55 km s _1 and compute the luminosity of pre- 
collision and post-collision stars as in Sect. [|. Then we 
add the luminosity of stars born in the central region and 
use the total luminosity in each annulus to calculate the 
resulting surface brightness, color, and index profiles. 
While the inner ring and the nucleus are currently HI 



poor with Shi < 0.3 Mq pc 2 (Higdon |l996|) , they have 



a significant amount of molecular hydrogen as sugg ested 
by recent CO observations of Horellou et al. ( 1998| ). We 
assume that molecular hydrogen in the Cartwheel is expo- 
nentially distributed in the central 9.5 kpc, with the peak 



surface density of 35 M© pc -2 in the nucleus and the scale 
length of 3 kpc (radius of the inner ring) . The total mass of 
H2 is then close 1.6 x 10 9 Mq, as estimated by Horrelou et 



al. (1998) for the solar metallicity. Metallicity in the cen- 
tral part of the Cartwheel is unknown. We assume that 
star- forming gas in the central 9.5 kpc is of solar metal- 
licity, as suggested by numerical modelling in Sect. |]. The 
star formation law in the nucleus and the inner ring is 
assumed to follow the classic Schmidt behavior: 



SFR 



a ! 



(6) 



with the rate of star formation a defined as in Eq. (|3j). 

We use two conventional models for star formation: 
a constant star formation and a short starburst. In both 
models star formation proceeds until 1.6 x 10 s Mq of stars 
is formed, which is 10% of the initial mass of gas. In the 
short starburst model we increase the rate of star forma- 
tion a to shorten the star formation event. 

Constant star formation in the nuclear region for the 
last 60 Myr. Assumption of the constant star formation 
activity in the Cartwheel's central region allows us to 
reach a better agreement between the model and observed 
B — V/V — K colors and Qbvk indices of the nucleus and 
the inner ring. The model Qbvk index of the inner ring 
lies within the observed limits, which is not the case in 
Fig. 5d. Now, however, the inner ring+nucleus and the 
outer ring have comparable Ha fluxes. This is problematic 
since observations indicate that the Cartwheel lacks strong 
Ha emission from the nucleus and the inner ring, while it 



has a large Ha luminosity in the outer ring (Higdon 1995, 
Amram et al. 1998). We assume that the lack of strong Ha 



emission from the nucleus and the inner ring is the result 
of a burst-like nature of star formation in the nuclear re- 
gion. If the active star formation in the nucleus-l- inner ring 
had ceased about ten million years ago, then the observed 
Ha flux would have fallen considerably at present. 

Short starburst in the nuclear region. Infall of gas into 
central parts of a galaxy is known to trigger starbursts. 
We assume that t he lar ge-scale radial infall of gas, re- 
ported by Higdon ( 1996 ), has triggered a short starburst 
in the Cartwheel's nucleus and the inner ring ten million 
years ago. The resulting model B — V/V — K radial color 
gradients and Qbvk radial profiles arc illustrated with the 
open circles in Figs. 7a and 7b respectively, while the filled 
triangles with error bars indicate the corresponding radial 
profiles observed in the Cartwheel galaxy. The sequence of 
Qbvk indices observed in the galaxy is well reproduced, 
as shown by the open circles in Fig. 7b. Now, the model 
peak value of Ha surface brightness at the outer ring is 
60 times larger than the model peak value at the nucleus, 
which is close to what is observed in the Cartwheel galaxy 
(Amram et al. |1998[ ). 

The difference of the model and Cartwheel's B — 
V/V — K colors of the inner ring (annulus I) and the 
nucleus, which is seen in Fig. 7a, can be explained in 
terms of internal extinction using the same approach as in 
Sect. ||. Implied by this difference, color excess E(B — V) 
of the Cartwheel's nucleus is m .43, which is equivalent to 
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E H i+h 2 of 19 M Q per 2 (Bohlin et al. |1978[ ). If molec- 
ular hydrogen is exponentially distributed in the cen- 
tral 9.5 kpc, then this value of £hi+h 2 might be below 
the detection upper limits of molecular hydrogen in the 
Cartwheel's nucleus. Hence, the internal extinction can 
be responsible for the difference between the model and 
observed B — V/V — K colors of the nucleus shown in 
Fig. 7a. This conclusion is also true for the Cartwheel's 
inner ring, where Ehi+h 2 = 4.5 M@ pc -2 is inferred using 
E(B — V) = m .l. Now, the young post-collision stellar 
populations and the old pre-collision stellar populations 
contribute equally to the total flux in the V- and iT-bands 
at the smaller radii (R < 3 kpc). 



6. Conclusions 

In this paper we model and analyse the B — V/V — K ra- 
dial color gradients observed in the Cartwheel ring galaxy. 
We use the Qbvk combined color indices to analyse the 
Cartwheel's B — V/V — K radial color gradients. Use of 
the Qbvk indices minimises the uncertainties in the ob- 
served B — V and V — K colors introduced by dust ex- 
tinction. We find that the Qbvk radial profile observed in 
the Cartwheel's disk falls naturally into three sub-profiles, 
which correspond to the nucleus+inner ring, the inter-ring 
region, and the outer ring of the galaxy. Combined color 
indices are insensitive to dust extinction and can be di- 
rectly understood in terms of stellar populations. This 
implies that the nucleus+inner ring, inter-ring region, and 
the outer ring of the Cartwheel might consist of different 
stellar populations. 

To model the optical and near-infrared color properties 
of the Cartwheel galaxy, we use a toy model, which imi- 
tates the propagating star formation in the galactic disk. 
Namely, we assume that the star formation rate is pro- 
portional to a Gaussian density enhancement, expanding 
radially in a gaseous disk. We use the population synthe- 
sis to calculate the color properties of stellar populations 
formed by this density enhancement. 

The results of color modelling strongly suggest that 
the pre-collision Cartwheel was a late-type spiral. The 
old stellar disk was embedded in an extensive gaseous 
disk of sub-critical surface density. Existence of the pre- 
collision stellar disk of low surface brightness fiy = 
23.5 mag arcsec -2 , extending out to the current location 
of the Cartwheel's outer ring at 16 kpc (Rq = 5 kpc), 
was first reported by Korchagin et al. ( 2001 ). However, 
we find that the pre-collision stellar disk might be smaller, 
extending out to « 12 kpc. Properties of the pre-collision 
stellar disk are typical for the late-type Freeman disks 
(Hy = 21.0 mag arcsec -2 and Rq = 2 kpc) rather than for 
low surfa ce brightness galaxies as found by Korchagin et 
al. ( 2001 ). The pre-collision gaseous disk had a mctallic- 
ity gradient ranging from z = z@/5 at the outer parts to 
z = zq in the central parts. Approximately 300 Myr ago, 
the Cartwheel galaxy collided face-on and near-centrally 
with a companion galaxy. This collision generated an ex- 



panding ring density wave, which triggered massive star 
formation along the wave's perimeter. 

The results of our modelling show that the young stel- 
lar populations, formed by the density wave, dominate in 
the outer ring. Thus, the wave of star formation has passed 
the extent of the pre-collision stellar disk and is currently 
moving in the predominantly gaseous, low-metallicity disk 
at the radius of 16 kpc. In the region between the inner 
and the outer rings, however, the young stellar popula- 
tions and the pre-collision stellar populations contribute 
comparable fluxes in JT-band. We find that the Cartwheel's 
B — V/V — K colors and Qbvk indices of this inter-ring re- 
gion cannot be successfully modelled without taking into 
consideration the relative input from the old, pre-collision 
stellar populations. 

The observed colors of the nucleus+inner ring are 
most difficult to model. Neither young stellar populations 
formed in the expanding density wave, nor their mixture 
with the old, pre-collision stellar populations can repro- 
duce the B — V and V — K colors of the Cartwheel's 
nucleus+inner ring. We find that an additional 10-Myr- 
ago burst of star formation in the nuclear regions, along 
with the visual extinction of Ay — l m .3, might be respon- 
sible for the peculiar colors of the Cartwheel's nucleus. 
This is not totally unexpected, as some observational ev- 
idence points to the existence of a low-level star forma- 
tion activity in the central regions (Amram et al. |1998 



Charmandaris et al. 199S ) and the large-scal e radi al infall 
of gas in the Cartwheel, reported by Higdon ( |1996| ), might 
have triggered bursts of star formation in the past. 
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